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Abstract 

A recently proposed statistical theory of the mean fields associated with the 
ground and excited collective states of a generic many-body system is ex- 
tended by increasing the dimensions of the P-space. In applying the new 
framework to nuclear matter, in addition to the mean field energies we obtain 
their fluctuations as well, together with the ones of the wavefunctions, in first 
order of the expansion in the complexity of the Q-space states. The physics 
described by the latter is assumed to be random. To extract numerical pre- 
dictions out of our scheme we develop a schematic version of the approach, 
which, while much simplified, yields results of significance on the size of the 
error affecting the mean fields, on the magnitude of the residual effective inter- 
action, on the ground state spectroscopic factor and on the mixing occurring 
between the vectors spanning the P-space. 
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I. INTRODUCTION 



We have recently addressed the problem of assessing the impact of randomness in the 
ground state of atomic nuclei In exploring this issue we assume random processes to 

be related to the strong repulsion the nucleons experience when they come close to each 
other inside the nucleus. This occurrence, while comparatively rare, is of crucial relevance 
for the nuclear dynamics and in the perturbative approach one tries to compute the related 
contribution, e. g., to the binding energy of the system, using G- matrix techniques. The 
latter, however, account solely for the two-body aspect of the short-range nucleon-nucleon 
(NN) correlations and do so only within certain approximations. 

Furthermore, the large violation of the Hughenoltz-Van Hove theorem occurring in the 
Brueckner-Hartree-Fock (BHF) nuclear matter calculations, — not to mention the well- 
known failure of the BHF theory in accounting for the empirical data, — clearly demonstrates 
that diagrams beyond the ladder ones are needed to satisfactorily describe the physics of 
atomic nuclei, in particular the short-range one. Higher order diagrams in the Brueckner- 
Bethe-Goldstone hole-line expansion have been computed, but they still appear to fall short 
of accounting for the experiments. 

Alternatively, one seeks for an exact solution of the nuclear many-body problem via the 
Green's function Monte Carlo method. This has actually been achieved lately in light nuclei, 
owing to the remarkable progress in the computational capabilities of modern computers (see, 
e. g., Ref. 1^. Yet, even the present days number crunching power is not enough to tackle 
the enormously complex problem represented by medium-heavy and heavy nuclei, precisely 
as it happens, in a different context, for the non-perturbative QCD. 

Moreover, given that these computational obstacles will eventually be surmounted, one 
will still have to face the question whether a non-relativistic potential model description is 
sufficient to encompass all the aspects of the nuclear many-body problem. Indeed, the issue 
of Lorentz covariance, of particular importance for the short-range physics in nuclei where 
large momenta are attained, is lurking in a corner, not to mention the ambiguities existing 
in defining and constructing a potential model description of the NN interaction at short 
distances, where it is so hard to disentangle particle and nuclear physics from each other 
and where off-shell effects become larger. 

It thus appears warranted to explore alternative routes in attacking the nuclear many- 
body problem. Hence, we have recently examined an approach to the theory of the ground 
state of atomic nuclei based on the concept of averaging rather than computing most of the 
physics related to the strong NN collisions at short distances. As a first step we have thus 
set up an energy averaging procedure suitable for constructing a mean field to be eventually 
identified with the shell model. Here, our procedure parallels the one successfully adopted 
in the development of the "optical model" of nuclear reactions to which it is in fact 
directly inspired. 

Clearly, for the mean field to have a meaning an assessment should be provided for the 
error associated with it. Indeed, if the error turns out to be large, then the concept of mean 
field is no longer tenable. Hence, we have derived an expression for the error through an 
expansion in the complexity of the states that the NN repulsion at short distances allows 
to reach. At the basis of the derivation lies the hypothesis alluded to in the beginning of 
this Introduction, namely that the matrix elements of the NN repulsion are random, thus 
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entailing the vanishing of the average value of the error, but, of course, not of its square. It 
is this one that our expansion (finite in finite nuclei and fastly convergent, as we shall prove) 
provides. 

The above outlined framework basically correspond to a statistical theory of the ground 
state (more generally, of the bound states) of the atomic nuclei. As a first example, we 
have attempted to implement it in the simplest among the latter, namely nuclear matter, 
assuming this system to represent, at least partially, the physics of heavy nuclei. Already 
in performing this task, far from trivial technical problems have been encountered, notably 
the one of performing the sum over the complex nuclear excited states building up the error 
(square of). We have been able to surmount this obstacle in Ref. 0. 

It should be furthermore observed that, at the present stage, even in nuclear matter our 
statistical theory cannot be made parameter-free. As a consequence, its predictive power is 
somewhat limited: Indeed, what we actually do amounts to correlate a set of observables of 
nuclear matter (or heavy nuclei), either measured or computed (via, e. g., the BHF theory), 
— like the binding and the excited states energies and the level density, — and "predict" 
on this basis the mean field energies of the ground and of the collective nuclear states, — 
together with the associated error, — the residual effective interaction, which, while well 
defined, is presently hardly computable, and in addition the ground state spectroscopic 
factor. 

Concerning the parameters entering in our approach, one is needed to characterize our 
energy averaging procedure, whereas the others are required whenever the experimental 
values of the above referred to observables are lacking or their theoretical evaluation not 
available. We shall discuss in particular the significance of the former, which plays a central 
role in our approach. 

This work is cast in the language introduced to build the unified theory of nuclear 
reactions [^, at the core of which lies the partition of the nuclear Hilbert space in the P 
and Q sectors. We implement this partition by inserting the random aspects of the physics 
of the nucleus in the latter. Hence the Q-space should not be viewed as specific of a given 
Hamiltonian, but rather it displays universal features and indeed it is in the Q-space that 
the energy average is performed. The P-space embodies instead the deterministic physics 
of the nucleus and its dimensions should be dovetailed for this purpose. This is why the 
chief scope of the present research is to expand the dimensions of the P sector of the nuclear 
Hilbert space, thus removing a major restriction of our past work, where the P-space was 
limited to one dimension only. 

Specifically, the enlargement of the dimensions of the P-space allows: i) to compute 
the fiuctuations of the nuclear wave function and not only of the energy, thus assessing the 
error associated both with the mean field energy and its wave function; ii) to show the fast 
convergence of the expansion for the error, essentially stemming from the rapid increase of 
the nuclear level density with the excitation energy; iii) to provide an improved estimate, 
with respect to our past work, of the mean field energy and of the spectroscopic factor of 
nuclear matter; iv) to predict a quenching of an order of magnitude of the matrix elements 
of the NN residual effective interaction as compared to those of the bare force. 

In the following we shall dwell on the above items shortly revisiting the theory in Sect. 0, 
deepening the meaning of the energy averaging in Sect. |T| and broadening the P-space in 
Sect. Next, we shall address the question of the spectroscopic factor and of the wave 
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function fluctuations in Sect. 0. Finally, the numerical results obtained on the basis of a 
schematic model will be presented in Sect. 



II. THEORETICAL FRAMEWORK 

Our statistical theory has been developed in Refs. flj-Q. For the convenience of the 
reader and to facilitate the understanding of the extension of the theory presented in this 
paper, we shortly recall in the following the basic equations characterizing our approach. 

As it is well-known, the splitting of the Hilbert space induced by the projection operators 
P and Q transforms the Schroedinger equation 

Hip = Etp (2.1) 

into the pair of coupled equations 

{E - Hpp){P^) = Hpq{Q^) (2.2a) 
iE-HQQ)iQiP)=HQpiPiP), (2.2b) 

the meaning of the symbols being obvious. From the above the equation obeyed solely by 
{Pip) is derived. It reads 

n{PiP)=E{PiP), (2.3) 

the P-space Hamiltonian being 

= Hpp + HpQ— — — 3^ Hqp, (2.4) 



i ) +Wqq 



with 



eQ = E-HQQ- Wqq (2.5a) 

Wqq = Hqp \ HpQ. (2.5b) 

rL/ — npp 

It is of significance that although Eq. ( p.3| ) is not an eigenvalue equation, since the energy E 
also appears in the denominator of its right hand side, yet its solutions only occur for those 
values of E which are eigenvalues of ( ^.1[ ) as well; in other words, a one-to-one correspondence 
between the values of E allowed by ( |2.1| ) and ( |2.3| ) exists (see next section for a further 
discussion of this point). 

Now, as already mentioned in the Introduction, in our framework the quantum deter- 
ministic aspect of nuclear dynamics is assumed to be embodied in the P-space, the chaotic 
one in the Q-space. Hence, the strategy of averaging over the latter follows, although, ad- 
mittedly, some fuzziness does affect this partitioning. To set up the averaging procedure we 
start by the recognition that the wave functions in the Q-space are rapidly varying functions 
of the energy E, viewed as a parameter classifying their ensemble. 

Accordingly, we average over this ensemble following the prescription 
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/oo _ 
dEp{E,Eo,e)f{E), (2.6) 
-oo 

/ being a generic function to be averaged over the variable E with the distribution 
p{E, Eq, e). The latter depends, beyond E, also upon the value Eq around which the average, 
— taken over a range of E essentially set by e, — is performed. A distribution convenient 
for our purposes is 

1 e*^'' 

p{E, Eo, e) = TT-Twp ^ ^2.7) 

zvrz E — [Eq ~ e) — 17] 

which is indeed correctly normalized being 

dEp{E,Eo,e) = 1 (2.8) 



(one should let 77 — 0"*" after the integration has been performed). Note that Eq. ( p.7| ) 
extends in some sense the Lorentz distribution of the optical model |^ to the situation of a 
zero width state. Hence the present formalism is especially suited to deal with the nuclear 
ground state, which is of course stable: We shall accordingly focus mainly on the latter in 
the following. 

We now recall that in [0-^ the Q-space wave functions was found to read 

(Q^) = —HQptfjo, (2.9) 
eg 

ipQ being an auxiliary function that in the end disappears from the formalism. By averaging 
Eq. ( |2.9| ) according to the prescriptions ( p.6| ) and ( |2.7| ), one then finds that the energy 
averaged wave function of the nuclear ground state in the P-space (here denoted by the 
angle brackets) obeys the equation 

H{Pij) = Eo{Ptlj). (2.10) 

In (|2.10|) Eq is the mean field energy and 

'H = Hpp + HpQ- —I Hqp (2-11) 

h) +Wqq{E = Eq) 



is the mean field Hamiltonian. This can be further elaborated since the singularities of the 
operator l/eg lie in the Imii^ < half-plane Accordingly, one gets 

(-) = —/ dE 



2m J-oo E-{Eo-e)-i7]E- Hqq - Wqq{E) 

^ Eo-e-HQQ-WQQ{E = Eo-e) Eo - e - Hqq - Wqq{E = Eq)' ^^'^^^ 

the last passage holding if the energy dependence of the operator Wqq is mild and if the 
parameter e is not too large (it should be not too small either, otherwise the energy averaging 
procedure becomes meaningless). 
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The insertion of ( |2.12|) into ( p.ll|) leads to the following useful alternative expression for 
the mean field Hamiltonian 



H = Hpp + VpqVqp 



1 



Eo-e-E' 



(2.13) 



where the energy dependent operators 



V, 



PQ 



H 



Eq — e — E 



\Eo-e-H, 



QQ 



(2.14a) 



and 



QP 



Eo-e-E 



\Eo-e-H, 



■H, 



QQ 



QP^ 



(2.14b) 



represent the residual effective NN interaction. The usefulness of the Eqs. ( p.l3| ) and ( |2.14| ) 
was realized in Ref . , where it was noticed that with their help the pair of equations ( p.2|) 
can be recast, as far as {Pip) is concerned, into the form 



{E-H){PiP) = VpQ{Qi,) 
{E-Hqq){Q,P) = Vqp{PiP), 



(2.15a) 
(2.15b) 



which is suitable for expressing the mean field fluctuations (the "error"). 

Indeed, by exploiting the spectral decomposition of the operator {E — H)^^ in terms of 
the eigenf unctions 0„ of the mean field Hamiltonian 7i, one gets from Eq. (p. 15a ) 



m = E 



E — En 

{<Po\vpQm) 



vpQm) 



1 



E-Eo \E-n 
which, upon left multiplication by (0o|, yields 

>o|Vpq|Q^) 



(2.16) 



(0O|P^) 



E-E^ 



(2.17) 



In the second term on the right hand side of Eq. ( |2.16| ), the prime stands for the omission 
of the n = term in the spectral decomposition. 

Next, the insertion of Eq. (|2.16|) into ( p.l5b|) leads to 



1 



where the operator 



E-hc 



Hqq + Vqp 



-VQp\^o){Mm, 



(2.18) 



E-n 



PQ 



(2.19) 
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has been introduced. Finally, by combining ( p.l7|) and ( |2.18| ), one arrives at the equation 



E — En 



o\VpQ 



E-k 



QQ 



-Vqp|0o), 



(2.20) 



which is the basis for computing the mean field energy error. 

Although Eq. ( p.20| ) is valid for any choice of the projectors P and Q, its use is in fact 
appropriate when the P-space is one-dimensional, as it was indeed the case in our past work, 
where this choice had been made for sake of simplicity. For example, for a two-dimensional 
P-space, as we shall later see, one should rather single out two, rather then one, terms in 
the spectral decomposition of the operator {E — 7i)~^ on the right hand side of Eq. ( |2.16| ). 

We refer the reader to Refs. for a discussion on how the average of the square of 
Eq. ( |2.20D is actually computed (the average of ( |2.2CI| ) of course should vanish) and on how 
the complexity expansion is organized. Here, we simply remind that in the present work, 
— as we have done in the past, — we shall confine ourselves to the leading term of this 
fast converging expansion. The difference between the present results and the previous ones 
thus stems entirely from the enlargement of the P-space. 



III. ENERGY AVERAGING 



To understand better the significance of the energy averaging distribution (|2.7|) let us 
study how it works in the simple cases of a bi-dimensional (A) and of a tri-dimensional (B) 
Hilbert space. 



A. Bi-dimensional Hilbert space 

Let Ixi) and IX2) be the two normalized states spanning the space. Here the only possible 
choice for the projectors clearly is 

and Q = |x2)(X2|- (3.1) 
Then, by expanding the operator 1/{E — Hqq), Eq. ( pl3|) can be recast as follows 



|P^) = 0, (3.2) 



which, upon multiplying from the left by and exploiting the idempotency of |x2)(X2|, 
simplifies to 



E - an- 



ai2 



{Xi\Pij) = 0, 



E — a22 

where the shorthand notations 

an = (Xil^^lXi), a22 = {X2\H\x2) and ai2 = {xi\H\x2) 



(3.3) 



(3.4) 
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have been introduced. This equation is trivially solved yielding the eigenvalues 



E4 



an + a22 ± V (an - ^22)^ + 4:\^i2 



(3.5) 



which coincide with those of H. It helps notice that the eigenvalues ( |3.5| ) are also found as 
intersections of the hyperbola 



E = aii + 



ai2 



uj — 022 



(3.6) 



with the straight line E = uj. 

Also the energy averaged Hamiltonian ( p.l3|) can be expressed in the basis spanned by 
Xi and X2 and one gets the mean field equation 



E - ail- 



ai2 



E — 022 

The latter is again trivially solved yielding 



(Xi|(P^)) = 0. 



(3.7) 



E. 



ail + 022 + e ± W (an - ct22 - e)^ + 46^ 



(3.8) 



which now corresponds to the intersections of the hyperbola ( |3.6|) (with E replacing E) with 
the new straight line E = oj + e. 

From Fig. |l|, where the solutions E± and E± are graphically displayed, it clearly appears 
that, while E_ = E_, the other solution E^ is much larger than E^, the more so the greater e 
is. It is thus clear that the averaging distribution (|2.7| ), while mildly affecting the eigenvalue 
of H lying in the P-space, drives away the one lying in the Q-space. 



B. Tri-dimensional Hilbert space 

The space is spanned by the normalized states IX2) and Ixa)- Now, two choices are 
possible for the projectors, namely 

P= |Xi)(xi| + |X2)(X2| (3.9a) 

Q = \X3){X3\ (3.9b) 

and 

P=\Xi){xi\ (3.10a) 

Q=\X2){X2\ + \X3){X3\- (3.10b) 

In both cases, Eq. ( p.3|) can be recast as follows 

{E - aii){E - a22){E - 033) - \ai2\^{E - 033) - \avi\^{E - 022) - \a23\^{E - an) = 0, 

(3.11) 
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(0 

FIG. 1. The eigenvalues of a bi-dimensional Hilbert space. The matrix elements of the 
Hamiltonian are taken to be an = —2, 022 = 2 and ai2 = 3, in arbitrary units. The exact 
eigenvalues (squares) and the ones of the energy averaged Hamiltonian 7i (circles) are shown. 
They correspond to the intersections with the straight line E = u: + e: The former with e = 0, the 
latter with e = 7. The stability of the lowest eigenvalue and the upward shift of the highest one 
are clearly apparent. 



which is the cubic equation yielding the exact eigenvalues. Note that Eq. (|3.11|) is easily 
obtained with the choice ( p.9|) , because in this case the operator {E — Hqq)^^ is expanded 
in terms of the idempotent operator |x3)(X3|- Not so with the choice ( p.lO|) , because now 



[E — Hqq)^^ should be expanded in terms of the operator ( p.lOb| ), which is not idempotent. 
Actually, the larger the powers of the latter are, the more cumbersome they become. Yet, 
also in this case it can be proved that Eq. ( |3.11D holds valid. 

Let us now examine the solutions of Eq. (|2.13|) : As in the previous bi-dimensional case 
it is convenient to display the solutions graphically. For the partition ( p. 101 ), one finds that 
they are given by the intersections of the e-independent curve 

^ , |ai2p(t^ — 033) |ai3p(u; — 022) , 0120*3023 , 0*2013023 /o io\ 

where 

V{u) = {u - 022)(CU - 033) - |a23|^ (3.13) 

with the straight line ^ = u; + e, as displayed in Fig. 0, where the case e = 0, — which 
clearly provides the exact eigenvalues Ei of the Schroedinger equation, — is also shown. 
From the figure, it transparently appears that Eq = Eq, whereas Ei» Ei and E2 » E2, 
the latter inequalities being stronger when the parameter e is large. 

In the case of the partition ( p.9| ), the solutions are given by the intersections of the curve 
obtained replacing 022 — O22 — e in Eqs. ( p.l2|) and (|3.13|) with the straight line E = lj + e, 



as displayed in Fig. ^. We face here a new situation, since now not only the straight line, but 
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K) 

FIG. 2. The eigenvalues of a tri-dimensional Hilbert space: The case of a one-dimensional 
P-space. The following matrix elements of the Hamiltonian are taken: an = —4, 022 = —2, 
033 = 3, ai2 = 2, ai3 = —5 and 023 = 3.5, in arbitrary units. The exact eigenvalues (squares) 
and the ones of the energy averaged Hamiltonian Ti (circles) are shown. They correspond to the 
intersections with the straight line E = lo + e: The former with e = 0, the latter with e = 10. 
The stability of the P-space eigenvalue and the upward shift of those belonging to the Q-space are 
clearly apparent. 




CO 

FIG. 3. The eigenvalues of a tri-dimensional Hilbert space: The case of a bi-dimensional 
P-space. The matrix elements of the Hamiltonian are taken as in Fig. ^. The exact eigenvalues 
(squares) correspond to the intersections of the straight line E = co with the continuous curve. The 
eigenvalues of the energy averaged Hamiltonian Ti correspond to the intersections of the straight 
line E = uj + e, with e = 10, with the dashed curves. Note the dependence upon e. The stability of 
the P-space eigenvalue and the upward shift of the one belonging to the Q-space are again clearly 
apparent. 
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also Eqs. (|3.12|) and ( p.l3|) are e-dependent. Yet, one again sees that for e = one recovers 
the eigenvalues Ei, whereas when e 7^ the intercepts occur for Eq = Eq and Ei = Ei, 
but for E2 » E2. Hence, we conclude that the action of the averaging distribution ( p.7| ) 
affects very little the eigenvalues belonging to the P-space, while pushing off the ones in the 
Q-space by an amount proportional to e. 



IV. ENLARGING THE P-SPACE 

As in our previous work, we address in the following the problem of the ground state 
of nuclear matter, again shortly summarizing below the equations derived in Refs. @J^). 
Indeed, the improvements brought in by the present research are more transparently appre- 
ciated by performing the analysis horning in mind the one formerly developed. 

As a preliminary in carrying out our task, it is necessary to define the operators P and 
Q. For this purpose the natural candidates as building blocks of the P operator appear to 
be the eigenstates of the mean field Hamiltonian ( p.l3|) , defined by the equation 



=E„|0„). (4.1) 

Their finding requires, however, the solution of a difficult self-consistency problem. Hence, 
in Refs. the simpler choice had been made of viewing as building blocks of P the 
Hartree-Fock (HF) variational solutions, which are trivial in nuclear matter. 

Furthermore, there, for sake of illustration, the severe limitation of a one-dimensional 
P-space was adopted by setting 

|Xhf)(xhf|, (4.2) 

Ixhf) being the HF ground state wave function of nuclear matter (the Fermi sphere). It was 
then shown that, on the basis of ( [4.2| ), one derives the mean field equation 



Eo = EnF + ^ , (4.3) 

Eq- e- E 

which relates the mean field (Eq), the HF (£'hf) and the true (E) energies per particle, and 
the fluctuation equation 




(4.4) 



where 



/3^= E l(^2p-2h|V^|XHF)P, (4.5) 
2p-2h 

the bras ('?/'2p-2h| representing the two-particle-two-holes (2p-2h) states of nuclear matter, 
whose average energy per particle is €2- We emphasize that all the quantities appearing in 
Eqs. (|4.3| ) and (|4.4|) are per particle, including the parameter e and the residual effective 
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interaction V, both of which should accordingly be thought of here as being divided by the 
nuclear mass number A. 

Note also that Eq. ( |4.4| ) gives the "error" of the mean field energy in the first order of 
the complexity expansion: Hence, it is contributed to only by the sector of the Q-space 
set up with the 2p-2h excitations. Actually, again for sake of simplicity, the states of the 
Q-space, although obeying well-defined, coupled differential equations (see Ref. 0), have 
been assumed to be adequately represented by the HF multi-particle-multi-hole solutions, 
an approximation not impairing the orthogonality constraint P ■ Q = 0. 

The sum in Eq. ([4.5[) is performed over the ensemble of the 2p-2h excited states lying 



in an appropriate energy range (in Refs. taken to be fixed essentially by the parameter 
e), whose number A/2 can be computed using the Ericson's formula for the density of the 
spin J A^-particle-A^-hole nuclear states, namely 

J) = Jluf.'u, . ^^L -P|-(2^ + ^f/ma% (4.6) 

where 



p\h\{N - 1)\ ^am^/^ 




9 = ^— and a2=^W--2, (4.7) 
2eF ^ ah 

T being the nuclear moment of inertia, ep the Fermi energy, a = A/S MeV~^ and E the 
excitation energy of the system. 

Now, Eqs. (|4.3| ) and ( [4.4|) , because of the double sign appearing in the latter, set up two 



systems, each one including two equations, in two unknowns. Two options are possible in 
selecting the latter: One can choose either the ground state mean field and true energies per 
particle, i. e. and — assuming the matrix elements of the residual effective interaction 
to be under control, — or Eq and when E is experimentally known, — which is indeed 
the case in nuclear matter. 

The latter is the course we followed (and continue to follow in the present work) requiring 
the coincidence of the two £^0 obtained by solving the two systems separately. Actually, and 
notably, both systems lead to the same formal expression for the mean field energy per 
particle, namely 

= i |(Ehf + £ + e) - v/(^HF-^-e)2 + 4/32| , (4.8) 

which holds valid for e < i?HF — E, the right hand side of this inequality being positive 
because of the variational principle. The above, when fP' —>■ 0, yields 

Eo = E + e, (4.9) 

in accord with ( |2.14| ), but also with ( [4.4| ), which, for V —>■ 0, gives Eq = E, entailing 
that when the residual effective interaction vanishes no fluctuations occur and hence the 
parameter e should vanish as well. 

On the other hand, the two systems yield two different expressions for the sum of the 
matrix elements of the residual effective interaction squared, i. e. 
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+1 



Ei{l + JN2/2) - (e + JN2/2I2) - EnF 



El -12 



Eiil + V-A/'2/2) - (e + JM2/2-t2) - EuF 



+ Ae{Ei - Euf) 



(4.10a) 



and 



'Eu — t2 



E„(-l + JM2m - (-e + JM2l1-t2) + E, 



mF 



+ 1 



E„(-l + VM/S) - (-e + VA/'2/2e2) + ^HF 



4e(E„-EHF)^ (4.10b) 



both of which vanish in the hmit e — 0, in accord with the previous discussion. 

Formula ( ^4.10a|) , solution of the first system of equations (the one with the "+" sign 
on the right hand side of ( [4.4D ), yields the value of the energy dependent quantity ( |4.5|) 
(here denoted by jSf) on the lower border of the energy band expressing the fluctuations of 
the ground state energy E (and thus encompassing ^o)- Formula (|4.10b| ), solution of the 
second system of equations (the one with the minus sign on the right hand side of ([4.4| )), 
provides instead ( [4.5|) (here denoted by /5^) on the upper border of the band (remember that 
E-e2< 0). 

Of course, of the energy per particle E we only know the experimental value, not the 
values on the borders of the band: As a consequence, we can only surmise the width W of 
the latter, thus providing two different inputs for the energy E appearing on the right hand 
side of ( |4.1U| ), namely Eu = E + W/2 and Ei = E — W/2. However, we can then explore 
whether, for a given W, a value for the parameter e can be found (not too large, not too 
small) such to have the two mean field energies per particle to coincide. If and when this 
search succeeds, then an orientation on W (or, equivalently, on the size of the fluctuations 
of the ground state energy) can be gained. 

It is of importance to realize that the very existence of this scheme stems from the 
occurrence of the same quantity jS"^ in both ( [4.3|) and ( [4.4|) . This, in turn, is due to an 
approximation introduced in Ref. It will be further discussed in Appendix ^ 

We now improve upon the above approach by first letting the projector P to encompass, 
beyond the ground state, the 2p-2h excitations of the HF variational scheme as well. Thus, 
instead of Eq. ([4.2|) , we shall write 



P = |xhf)(xhf| + IXhf)(Xhf 



2/3 



(4.11) 



where the sum is meant to be extended to the whole set of 2p-2h HF excitations |xhf)- 

With the choice ( [4.11| ) the mean field Hamiltonian ( p.l3| ) is then defined and one can 
compute the mean field ground state energy per particle, Eq = (0o|'^|0o), using for the ket 
|0o) the expression 



d) = SoIxhf) + 5I*2|Xhf) 



(4.12) 
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and accounting for the influence of the Q-space on the ground state mean energy per particle 
in first order of the complexity expansion, i. e. by setting 



Q = EIXh^f)(Xh'^fI, (4.13) 



where the sum runs over the whole set of the HF 4p-4h excitations. In (|4.12|) sq and S2 



are complex coefficients, fixed, in principle, by Eq. ( [4.1| ) and satisfying the normalization 
condition 

koP + El4r = l- (4-14) 



After straightforward, but lengthy, algebra using ( [4.12| ) and ( [4.13|) one finally arrives at the 
following new mean field equation 

Eo = \so\'EnF + E l4r(XH"Fl^lxl"F) + 2.^ E 4{XnF\V\x'^F) + 

7 7 

+ J, ^ j, j:\4\'\{xt^,\V\xl^F)\\ (4.15) 



V being the bare NN potential. 

Notably, the above equation turns out to formally coincide with ([4.3|) . Indeed, one 
recognizes in the first three terms on the right hand side of Eq. ( |4.15| ) just the mean value 



of the original bare Hamiltonian H in the state |0o); given in (|4.12|) . In the thermodynamic 
limit characterizing nuclear matter what survives out of these pieces is just the HF energy, 
since the correction to the latter due to a finite number of particle-hole excitations goes to 
zero when divided, — in order to get the energy per particle, — by the infinite number of 
particles in the system: In other words, the 2p-2h admixture into ([4.12|) does not change the 
expectation value of the Hamiltonian. 
Hence, defining 

e = Y.W2?\{xZ\V\x'^^)\\ (4.16) 

Pi 

( [4.15| ) can be recast into the form 

E, = E^^+ ^] (4.17) 

whose similarity with Eq. ( |4.3|) is transparent. 

The same will take place for any admixture of A^p-A^h HF excited states in |xhf): Hence, 
in our framework different choices of the projection operator P lead to the same structure 
for the mean field equation for nuclear matter. This invariance does not hold in finite nuclei. 

The only, of course important, difference between (13) and ( 4.16 ) is that in the latter 



the residual interaction V induces transitions from 2p-2h to 4p-4h states, rather than from 
the Fermi sphere to the 2p-2h states. 

Concerning the fluctuation equation one can again use (|2.20|) , with the state |0o) given 
now by Eq. ( |4.12| ). Then invoking the randomness of the phases of the Q-space wave 



14 



functions (RPA) and proceeding exactly as previously done in one deduces the new 

fluctuation equation 




E-Eo = ±-^J^C, (4.18) 



where denotes the average energy per particle of the 4p-4h HF states. In (|4.18| ), A/4 



represents the number of 4p-4h excitations contributing to the sum over the index /3 in 

(EM)- 



Hence, the "formal invariance", with respect to the choice for the projector P, holds for 
both the equations at the core of our statistical approach to nuclear matter in first order of 
the complexity expansion. 

Indeed, the inclusion of A^p-A^h states (with N > 2) into the P-space would merely imply 
the replacement, in ( |4.18|) , of Af^ with Nn+2 and, at the same time, to have defined in 



terms of the matrix elements of V between iVp-iVh and (iV+2)p-(A^+2)h states. In addition, 
one should of course insert in the energy denominator the average energy per particle of the 
{N + 2)^-{N + 2)h HF states. 

Therefore, in this connection one clearly sees that the extension of the P-space rapidly 
leads to the vanishing of the fluctuations, owing to the very fast increase of the number Mn- 

However, it should be realized that Eq. ( |2.2CI| ), as it stands, is not in general a good 



starting point to derive the fluctuation equation. Indeed, it selects out only one term in the 
spectral decomposition of the operator 1/{E — H), which, as already pointed out, is justified 
for a one-dimensional P-space, but not so for a multidimensional one. Hence, we depart 
from our past work [0,^, not only by employing Eq. ( [4. Ill ), but also by writing, in place of 



(El). 



IP,) . ^ 1^.) <M!^ ^ ( 1 )' (4.19, 

11/ — ii/Q III — 1I12 Kill — 11 J 

with an obvious meaning of the double primed operator in the last term on the right hand 
side. Clearly, Eq. ( [4.19|) does not follows directly from ( |4.11|) , but it assumes that in the 



spectral decomposition of the {E — Ti.)^^ operator only one prominent (collective) state |02) 
enters beyond the ground state |0o)- 

Sticking to this model, instead of Eq. (|2.18|) , we likewise use 

IW = ^ \^2) VQp{\<Po){Mm + |02)(02|P^)}, (4.20) 

being 

h^QQ = Hqq + Vqp (-^]" VpQ. (4.21) 



.E-n 

Then, by left multiplying ( [4.19| ) (with E = Eq) by (0o| and using ( [4.20| ), we obtain 



-^0 — Eq — hQQ \<Po\-t^Wi 
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which generahzes Eq. ( p.20|) . 

In a perfectly analogous fashion, by left multiplying ( [4 .191) (with E = E2) by (02 1 and 
using ( ^.201 ), we obtain 

In the above, Eq and E2 stand for the first two exact eigenvalues of the Schroedinger equation; 
Eq and E2 for the corresponding quantities associated with Eq. ([4.1D. Of course, Eq. ( |4.23| ) 
has no counterpart in our previous work: Here, it shows that, in our frame, all the energies 
of the P-space fluctuate. 

Now, the energy averaging of ( [4.22| ) and ( [4.23| ) vanishes by definition, but the energy 
averaging of their square, which yields the "error", does not. Hence, proceeding along the 
lines of Refs. , we subtract on the right hand side of both equations their average values, 
square the expressions thus obtained and make use of RPA, keeping of our expansion in the 
complexity of the Q-space states the first term only. Next, we exploit the structure of \4>2), 
which, like \4>o), must be normalized, orthogonal to |0o) and of the form 

l02) = EcflXH^F) + co|XHF). (4.24) 
We thus finally arrive at the equation 

Eo-Eq = ±Tr^\[^iC" + ^e^)' (4-25) 



where the further definition 



e = T.^Vcl\{xZ\V\x^^)? (4.26) 

/37 



has been introduced in addition to ( [4.16|) ; moreover, we have set 



(4.27) 



(0o|P^)' 

Likewise, for the energy of the 2p-2h state of the P-space one obtains the fiuctuation equation 



where, naturally. 



E2-E2 = ±-^J-^{r^' + ^), (4.28) 
E2 - ei \ Mi r 



v' = Y.\cl?\{xZ\V\x%)?. (4.29) 

/97 



We thus see that the fiuctuation equations ( |4.25D and ( [4. 281 ) are actually coupled through 
the term ( ^1261) . 
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Concerning the mean field equations, it is clear that, in the enlarged scheme brought 
about by the projector ( 4.11| ) (actually, by our drastic modelling of the solutions of the 



mean field Hamiltonian ( p.l3| ) with P given by ( [4.1 1| )) and by the \Pil)) given in ( |4.19| ), an 
equation should exist also for the energy of the 2p-2h state. It can be derived by computing 
E2 = {<t>2\'H\4'2) and, notably, it turns out to read 

E2 = + ^ f p , (4.30) 

(2) 

-^HF representing the HF energy per particle of the system in the 2p-2h excited state. To 

(2) 

identify this one is not necessary, since we split E^^p into a part associated with the HF 
ground state and a part associated with the 2p-2h excitation energies, both per particle. 
Since the latter vanishes in the thermodynamic limit, — as previously noted in commenting 
Eq. ( [4.15| ), then E^l = Eb^f , a relation we expect to be approximately fulfilled also in a 



heavy nucleus. 

We conclude from the above analysis that our approach leads to a set of mean field 
equations, one for each of the states lying in the P-space: These equations, unlike the 
fluctuation ones, are not coupled. 

V. NORMALIZATION AND FLUCTUATION OF THE P-SPACE GROUND 

STATE WAVE FUNCTION 

In our framework one could identify the ground state spectroscopic factor 5* with the 
square root of the norm of {Pip), the system's ground state wave function projection in P- 
space. To work out an explicit expression for S we commence by exploiting the completeness 
of the normalized eigenstates of H. Hence we write 

M 

^ {P^IP^P) = 5:(P^|0„)(0„|PV^) = 1 - {QijlQij). (5.1) 

n=0 



Now, if we confine ourselves to set M = 1, as in Refs. then the Eq. ( p.l8| ) for \Qip) is 
warranted and we rewrite ( p.l| ) as follows 

S' = l~ (0o|V-pQ \ VQp|0o)|(0o|PV^)r. (5.2) 

Moreover, since the choice M = 1 clearly entails 

S' = \{<Pom\^ (5.3) 
by exploiting ( |2.20|) Eq. (|5^ ) can be recast into the form 



dEo Eq — e — Eq 

Finally, employing ([4.8D the expression, previously obtained in 



(5.4) 
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3 lEuF- Eo-e- 2dp^/dEo e 

2 2 . /(Ehf - ^0 - e)' + 4/52 Eo-e-Eo 



(5.5) 



follows, in which the energy derivative of the sum of the square moduls of the vacuum-2p- 
2h matrix elements of the effective interaction appears (its explicit expression is given in 
Ref. [Q). Note that ( |5.5| ) goes to one as /J^ —>■ 0, as it should. 

If, however, the expression ( [4.20[ ) for the Q-space wave function is used, then, rather 
than (|5.2|) , one obtains for the ground state spectroscopic factor 



-{<Po\Pi^){Hm*{HVpQ ^ 



Vqp\4>o)-\{(P2\P^)\'{4>2\Vi 



PQ 



(Eo - h'ahY 



jiy-yQp\<p2) 



(2) N2 



^QP|0O) 



{Eq - hgq) 

[Eq - iIqq) 

The above, with the help of Eqs. ( f4.22| ) and ([4.23|) , and assuming 

1 _ 1 

Eq - h^QQ E2 - h^Q^Q 

can be further elaborated and one ends up with the expression 



(5.6) 



(5.7) 



1 - 



d{E2-Eo) E^-Eg 
dEo ^ Eo~e~ Eo 



\{h\Pi^)? 



1 + 



dE. 



(5.8) 



which reduces to ( |5.4|) , as it should, if {(j)2\P'ip) — ^ 0, i. e. for a one- dimensional P-space. 
Since, from Eq. (|4.27|) , 



l(02|P^)P 



1 + r2 



then Eq. ( ^.8|) can be recast as follows 



1 + 



1 



Eq — e — Eq 



e + 



1 + 



(^E2 — Eq 



1 (dEQ ^dE, 
l + r2 \dEQ dEQ 



(5.9) 



(5.10) 



which again reduces to (|5.4| ) as r — > 0. 

We now address the problem of the fluctuation of \Pip). For this scope, we focus on the 
ground state and, by combining Eqs. (|2.10| ) and (|2.15|) , we obtain 



{E - nw^) - \{pm + {E- Eq)\{p^)) = VpqIq^) 

(the angle brackets meaning energy averaging). 



(5.11) 
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Then, if use is made of the expression (|2.18|) for and of the spectral decomposition 
of the operator {E — H)'^, one gets 

+-^\<Po){<Po\VpQ—^-—VQp\<j),){<j)o\P^) 

^ — ^0 ^ — iT-QQ 

-{E - Eo) E ^r^l0n)(0n|(P^)). (5.12) 

Now, since \{Pil))) oc |0o), from the above finally it follows 

1 \'. 1 



|p^) - \{Pi,)) = (1 - i0o)(0oi)-' (e^) ^^o:^37^^q^I'^o)(<^oI^^) 

= (0o|P^) (^-^)VpQ^^-— Vqp|0o), (5.13) 



which vanishes when P is given by (\i.2[ ), since clearly |PV) does not fluctuate in a one- 
dimensional P-space. 

If, on the other hand, P is given by Eq. ( ^.11| ), then the above can be computed, in first 
order of the complexity expansion, using ( [4.12| ) for |(/)o), ( ^4.24|) for |02) and ( [4.13| ) for Q. 
One ends up with the expression (we set E = Eq to conform to previous notations) 

[I^V^) - liPMi = {MP^)j}^ E {x^.\v\xt^.)^^{xt^F\v\xZ), (5.14) 

^0 - ^2 i3f^,^ Eq - e^p 

which can be further simplified invoking the randomness of the phases of the wave functions 
in the Q-space and again introducing the average energy per particle q for the 4p-4h HF 
excited states. Then, with the help of Eq. ([4.26|) , the formula 

is derived. It gives the fiuctuations of the wave function in first order of the complexity 
expansion. 

Notice that here the scalar product {(f)o\Pip), unlike in Refs. does not coincide with 
the spectroscopic factor S, as defined in (lOp , since Eq. ( |5.15|) has been obtained using an 
enlarged P-space. Rather, it measures the amount of the true ground state wave function 
of the system embodied in the mean field state |0o)- 



VI. RESULTS FROM A SCHEMATIC MODEL 

The numerical implementation of the formalism described in the previous Sections rep- 
resents a challenging task. 

Therefore, as a first step toward this goal, we develop here a version of our approach 
that, while simpler, has the merit of being amenable to numerical results. 
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To start with, because to solve Eq. ( |4.1| ), with the projector P given by ( [4.11|) , is a major 
task, we assume the two lowest solutions of the latter to be of the form 

I0o) = SoIxhf) + S2\Xrf^^) (6- la) 

and 

102) = -S2|xhf) + SoIxbf^)^ (6-lb) 

where Sg + S2 = 1. Of course, S2 remains to be fixed. 

This assumption corresponds to take in the sum on the right hand side of ( 4.11| ) one 
term only and places all the strength of the 2p-2h eigenstates of in a single mode. Then, 
by defining 

^^ = El(XH'^Fl^lxg'/')r, (6.2) 
7 

the expressions ( [4.16| ), ( |4.26| ) and ( |4.29| ) can be recast as follows 

= slfi'^ (6.3a) 



e = S2^1 - slfi^ (6.3b) 

= (1 - sDfi^ (6.3c) 

where, beyond /i^, only the coefficient S2 enters. 

Accordingly, the mean field and the fluctuations equations now read 

2 2 

Eo = EuF + ^ ^ (6.4a) 

E2 = E^^ + ^ — (6.4b) 



Eo — e — En 



and 



Eo-Eo = ^xi^^-^^i'l ± ^^2^/1 - si) (6.5a) 




E2-E2 = ±J-^-^-^(l ~sl± -S2^l - si), (6.5b) 



where E^^f and 64 have been previously defined. We thus see that each dimension added to 
the P-space entails the occurrence of two new systems in two equations. Furthermore, the 
double sign inside the round brackets on the right hand side of Eqs. ( |6.5|) corresponds to the 
two options one has in choosing sq, namely 



So = ±Vl-si. (6.6) 

Although the above equations are strictly valid in the thermodynamic limit only {A — >■ 
00), we shall make an heuristic use of them also for the nucleus A = 208, which is of course 
finite, but large. 
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Concerning the ground state spectroscopic factor ( ^.10|) , in the present model the energy 



derivative appearing in S"^ is easily found (from ( |6.4|) ) to read 

dEo^l^ l EnF-Eo-e-2d{slfi^)/dEo 

where account is taken of the energy dependence of the coefficient S2. 

In accord with the general theory previously discussed, both Eq and E2 lie inside an 
energy band expressing the "error" they are affected by. In conformity, also the ground 
state spectroscopic factor should be computed on the lower and upper borders of the ground 
state energy band: Hence, we need to compute on both boundaries Eq. ( |6.7| ), which in turn 
requires the knowledge of the energy derivative of the quantity (s2/i)^. 

For this purpose, in principle one should use the residual effective interaction on the 
lower boundary, which is found to be 

4f^l = do \ ^ {^0(1 + do)-t- dolA - EnF 

+ \l[Ei{l + rfo) - e - dot^ - Ehf]2 + 4e(E^ - Ehf)}, (6.8a) 
and on the upper one, where it reads 

4^1 = do °^ ^ {Eo{-1 + do) + e- d^e^ + ^hf 

+ V'[Eo"(-l + do) + e - doe^ + EnF? + 4e(Eo" - Erf)}, (6.8b) 



having set 




do = \h^ , (6.9) 



(the double sign stemming from (|6.6|) ), and where, as usual, Eq and Eq are the true ground 
state energies on the lower and upper borders of the band: As already discussed, these are 
expressed in terms of the band's width W . 

However, in practice, owing to the energy dependence of both the residual effective in- 
teraction /i and the coefficient S2 appearing in the ground state mean field wave function, 
the energy derivative of ( |6.8|) turns out to be very cumbersome to carry out. Accordingly, 
we evaluated numerically all the derivatives with respect to the exact ground state energy 
£^0 entering into the expression ( |5.1CI| ) for S'^. 



Note that d^ —>■ yA/'4/2 when r — »• 0: One thus formally recovers in this limit the results 
of a one-dimensional P-space [^. 

In solving the systems of equations ( |6.4a{ )-( |6.5a ) and ( |6.4b| )-( [6.5b| ) we have chosen in 



( |6.6|) the plus sign. The other option need not be considered because it will provide the 
same results, but for the change of sign of S2. Indeed, our equations are clearly invariant 
under the simultaneous transformations Sq —Sq and S2 — >■ — S2, as it should be as the 
latter simply correspond to a change of phase in the states (|6.1|). 
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Concerning the other quantities needed in our equations, we have assumed for the binding 
energy per particle of nuclear matter 



Eo = -16MeV (6.10a) 

for the ground state and 

E2 = -14.8MeV (6.10b) 

for the 2p-2h "collective" excited state. The latter value is hinted by the calculations of 
the (7,nn) and (7,pp) photo-absorption cross-sections in nuclei (n being a neutron, p a 
proton). 

For the HF energies per particle we have taken 

EHF = -12MeV (6.11a) 
= -11.96 MeV (6.11b) 
Q = -11.92 MeV. (6.11c) 

The first of the above actually corresponds to a typical BHF outcome (see, e. g., Ref. [0]). 
The other two stem from the minimal energy (= 4 MeV) of a ph excitation in a Wood-Saxon 
well with parameters adjusted for the nucleus Pfe^os 
The choice of a minimal ph energy has been made: 

a) because the only solution of our schematic model obtains with the residual effective 
interaction weakened by an order of magnitude with respect to the bare one: Hence 
the latter can only effectively connect with the lowest 4p-4h excitations (see the results 
below); 



b) because only with the minimal excitation energy of the 4p-4h states Eq. ( f4.6|) provides a 
number of states in an energy range set by e still compatible with non-trivial fluctuation 
equations (|6.5|). 



In connection with the point b) we have limited ourselves to use formula (|4.6| ) for 4p- 
4h excited states with angular momentum J = 0. It should also be added that one can 
legitimately question the validity of (4^), which is deduced in the framework of a Fermi gas 



by purely combinatorial methods. 

With the energies per particle above quoted, — which should clearly be viewed as merely 
orientative but for ( |6.10a] ), — we have solved the pair of systems corresponding to Eqs. ( |6.4aj ) 
and (|6.5a|) and, as well, the one corresponding to Eqs. ( |6.4b| ) and ( |6.5b|) , fulfilling the equal 
mean field constraint in both cases. In addition, for sake of simplicity, we have required the 
nucleons' residual effective interaction yU to come out the same in both instances in accord 
with the schematic nature of our model. 

Our findings for Eq, E2, fif and /x^ and for the ground state spectroscopic factor are 
displayed in Table |. These results have been obtained for a specific choice of e, which 
fulfills the inequalities < e < -Ehf — E = 4 MeV, S2 and r. Notably, in order to find a 
solution one has to restrict these parameters to a rather narrow range around the values 
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Eo (MeV) E2 (MeV) 






S2 e (MeV) r 


-16.001 -15.901 


25 26 


1.016 1.000 


0.71 3.2 1.02 



TABLE I. The solution of the schematic model 



quoted in the Table. Also the band width W cannot be much larger then 0.1 MeV in order 
to get a solution. The numbers in Table | correspond to the choice W = 0.1 MeV. 

Although we consider the outcomes of this model as merely orientative and without any 
pretense of being in touch with the real physics, both because of the crudeness of the model 
itself and because of the uncertainties affecting our inputs, yet they are worth a few remarks. 

First, the eigenstates of the mean field Hamiltonian 7i appear to have equal projections 
on the two vectors spanning the bi-dimensional P-space of our model, as it follows from the 
values we have found for S2. 

Second, it is clear that, with respect to our past work with a one- dimensional P-space, a 
reduction of the error affecting the mean field energy was to be expected. It turned out to 
amount to an order of magnitude. This fact points to a fast convergence of the complexity 
expansion or, alternatively, to a contribution of the 6p-6h excitations to the binding energy 
per particle (in nuclear matter) of less then 1%. This estimate is of course very rough: In 
fact, the magnitude of the error arises not anly from the dimensions of the P-space, but 
as well from the proximity of the HF (or BHF) solution to the experimental value of the 
binding energy per particle. 

Third, the massive reduction of the residual effective interaction as compared to the 
bare one, already found in Refs. is here confirmed. Unfortunately, in our past work 

we had overlooked a factor of A (the interaction under consideration is per particle) : The 
quenching there predicted was accordingly wrong, but still should be quantified in one order 
of magnitude (see Appendix |BD. 

Finally, a comment is in order on the spectroscopic factor S. For this purpose we write, in 
terms of the states |xhf) and Ixhe^*^) (0^5 alternatively, |0o) and |{/)2)) spanning our P-space, 

\Pi^) = IXhf)(xhf|V') + IXhf^^)(Xhf^''|V') 

= |0O)(0O|PV^) + |02)(02|P^). (6.12) 

From the above and (|6.1| ) it follows easily 

(6.13) 

{(j)2\P^) = -{XUf\iP)s2 + (XhF^^IV^)so- 

(6.14) 

Inverting these equations and exploiting ([4.27|) , we obtain 

(XhfIV^) = (001^^^)50 - {<P2\P^)S2 = (0o|P^)(so - rs2) (6.15a) 
(Xhf'V) = {MPi')s2 + {(p2\P^)so = {MP^)is2 + rso). (6.15b) 

Hence, with the values for S2 and r quoted in Table |, we conclude that 
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(XhfIV^) = 0, (xhf'V) = 1, (6.16) 

an outcome we are hardly ready to view as realistic. 

On the other hand, taking the minus sign inside the round brackets in the fluctuation 
equations ( |6.4| ) and also in ( |6.9| ) (which amounts to the substitution S2 —S2), we get 

(XhfI^) = 1, (Xhf'V) = 0, (6.17) 

which, while drastic , appears to provide a physical picture more consistent with the simple 
model we are employing. In fact, the above results hold because {(polPip) = 1, which follows 
from the normalization condition 

K0o|P^)r+|(02|P^)|' = l (6.18) 



and from the finding that r = 1. Clearly, both ( |6.16D and ( |6.17D lead to a value of one for 



the spectroscopic factor S, that is the result we have found. 

The two radically different solutions (|6.16|) and (|6.17|) can be reconciled through 



Eq. ( p.l5| ). The latter, indeed, when implemented in our schematic model, yields for the 



fluctuations of the ground state wave function 

[|^^)-|0o)]i = |02). (6.19) 

Likewise, by either simmetry arguments or explicit calculation, one can get for the collective 
2p-2h excited state 

[|P^2p-2h)-|02)]l=|0o). (6.20) 

Thus, it appears that, in the schematic model presented in this Section, while the fluctuation 
of the energy is puny, the fluctuation of the wave function is the largest possible. This 
correlation is worth being explored in more realistic models. 

In conclusion, it appears warranted to view the scalar product {(pQlPip) — which is 
informative on how much of the true ground state wave function of the system is embodied 
in the corresponding mean field one, — as a better representative of S: In this last instance, 
we would obtain S = 0.7. 



VII. CONCLUSIONS 

In this paper we have first shortly revisited the statistical theory of the ground and 
excited collective states of general many-body systems developed in Refs. DTj-Q. 

Next we have presented an improved version of the approach by allowing for a P-space 
with larger dimensions than in Refs. [0,0. In applying this new framework to nuclear 
matter we have obtained both the mean field energies and the fluctuations (the "error" ) of 
both energies and wavefunctions. This has been done in first order of the expansion in the 
complexity of the states of the Q-space, which, in our scheme, hosts the random aspects of 
the nuclear dynamics. Notably, the extension of the P-space allows to acknowledge the fast 
rate of convergence of the complexity expansion, due to the rapid growth of the nuclear level 
density with the excitation energy. 
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Finally, as a first step toward the testing of the theory on the physics of the atomic nuclei 
we have worked out a schematic model that renders the approach amenable to numerical 
predictions. While rough, the model has nevertheless helped us in clarifying basic aspects 
of our statistical approach. 

In conclusion, we feel the duty to state that the statistical approach here presented ger- 
minated out of the ideas discussed by Herman Feshbach in Ref . |l|] . It has been subsequently 
pursued by him in collaboration with us in Refs. |^,^ and in this work. 
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APPENDIX A: 

The derivation of the coupled equations at the basis of the present approach is possible 
because the same quantity (3"^ occurs in both ([4.3|) and ( [4.4|) . This occurrence stems from 
an approximation introduced in Ref. H], namely 



E((^2/3|V^|0O)*(^2,|V'|0O))((^2/?|^|<^O)(^2,|^|</'O)*) 



^2A\Y, ((^2/3|V^|0o)*(V^2,|V^|0o)) E ((^2;3|V^|0o)*(^2,|V^|0o)). (A1) 
1/37 J /37 

In the above, the operator A corresponds to the arithmetic mean of the quantity within 
the angle brackets over the set of the 2p-2h states in the energy interval e centered around 
e2 (actually, one should take the real part of the right hand side of the previous equation: 
However, for the sake of the argument we shall consider real quantities). 

This approximation gives reasonable results if the matrix elements are roughly of the 
same order of magnitude: This is the reason why the sum in (|A1D is limited to an energy 



interval around the mean 2p-2h energy where the matrix elements are relevant. 
We give here a few calculable examples that prove the above statement 

iVEliSinH^ 



2Et,sinH2]2 



N 



iVEliexp(-n) 



0.725, 


= 


10 


0.746, 


= 


100 


0.750, 


N = 


oo 


0.636, 


N = 


10 


0.663, 


N = 


100 


2/3, 


N = 


oo 


2.311, 


N = 


10 


23.11, 


N = 


100 


oo. 


N = 


oo 



2[Eliexp(-n)]2 - 
Were the approximation exact, then the above ratios should evaluate to one. 



(A2) 

(A3) 
(A4) 
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APPENDIX B: 



We report in the following the expression for — the sum of the square of the vacuum 
— * 2p-2h matrix elements, — for the schematic interaction employed in Refs. [0,0 , that is 

Vir)=gA- 9B- (Bl) 

r r 2 

Px being the exchange operator in coordinate space. 

Following the Appendix B of Ref . , can be recast as the combination of direct and 
exchange contributions such as 



,2 Q \ 1 

ir — ^ 

9 17, _o ~o\ l 



= ^77— <{ — ifJ'A arctan/iA - /^s arctan/xs) (B2) 



4:kF [flA-i^B 



1 1 

+ 



/i^(3 + arctan ( ) — fJ'A^'^ + /^a)^ arctan ( 

Ki^bJ KfJ^A, 



(B3) 



and 



Pl^^ = A—^J dq _^ q2 J dr—^sm{qr)^[sm{kFr) — kFrcos{kFr)]'^ 

+ [sin(A;i7r) — sm{r^Jkp — q^) — kpr cos(fci?r) + r \J kp — q^ cos(r-y/ k'jp — g^)]^!, (B4) 

where jlA,B = f^A,B/'^kF and g = 2kp/37i^. 

In the latter, a minor error in the expression reported in Ref. 0, essentially of no 
numerical consequence, has been corrected. More importantly, the necessity of normalizing 

to the mass number A had been overlooked in Refs. 0,0: Hence, the previously found 
quenching of the residual effective interaction with respect to the bare one should also be 
reduced by a factor A. For example, taking ^osp]^ orientation, the quenching would 

turn out to be of an order of magnitude. 
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